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ABSTRACT 

Motivated by the measurement of redshifted Fe Ha lines during type I X-ray bursts on EXO 0748-676 
(Cottam, Paerels & Mendez), we study the formation of atomic Fe lines above the photosphere of 
a bursting neutron star (^T e g w 1 — 2keV). We discuss the effects of Stark broadening, resonant 
scattering and NLTE (level population) on the formation of hydrogenic Fe Ha, Lya and Pa lines. From 
the observed equivalent width of the Fe Ha line, we find an implied Fe column of 1 — 3 x 10 20 cm~ 2 , 
which is 3-10 times larger than the Fe column calculated from the accretion/spallation model of 
Bildsten, Chang & Paerels. We also estimate that the implied Fe column is about a factor of 2-3 
larger than a uniform solar metallicity atmosphere. We discuss the effects of rotational broadening 
and find that the rotation rate of EXO 0748-676 must be slow, as confirmed by the recent measurement 
of a 45 Hz burst oscillation by Villarreal & Strohmayer. We also show that the Fe Lya EW » 15-20 
eV (redshifted 11-15 eV) and the Pa EW w 4-7 eV (redshifted 3-5 eV) when the Ha EW is 10 eV 
(redshifted 8 eV). The Lya line is rotationally broadened to a depth of w 10%, making it difficult 
to observe with Chandra. We also show that radiative levitation can likely support the Fe column 
needed to explain the line. 

Subject headings: stars: abundances, surface - stars: neutron - X-rays: binaries, bursts - lines: 
formation 



1. INTRODUCTION 

The recent observation of gravitationally redshifted 
atomic transition lines during type I X-ray bursts on the 
low-mass X-ray binary, EXO 0748-676 (Cottam, Paerels 
& Mendez 2002, hereafter CPM) opens a new window 
on neutron star (NS) structure. Using 335 ks of XMM- 
Newton observations during its calibration and commis- 
sioning phase, CPM collected 3.2 ks of Type I bursts 
which was split into a high temperature part (fcsTBB ~ 
1.8 keV) and a low temperature part (IcbTbb < 1.5 keV). 
In the high temperature part they identified the Ha 
(n = 2 — > 3) transition of hydrogenic Fe with a red- 
shift of z=0.35. The same transition was identified for 
He-like Fc with the same redshift in the low temperature 
part. Though this is not the first claim of an absorp- 
tion feature during type I X-ray bursts, previous Tenma 
and EXOSAT observations of a 4.1 keV absorption fea- 
ture on 2S 1636-536 (Waki et al. 1984; Turner & Bree- 
don 1984), EXO 1747-214 (Magnier et al. 1989) and 
X1608-52 (Nakamura, Inoue & Tanaka 1988) were never 
confirmed with more sensitive instruments. CPM's ob- 
servation is unique in that two lines in distinct states of 
the type I burst with the same redshift were observed, 
implying that we are seeing the NS surface. 

In our previous work (Bildsten, Chang & Paerels 2003, 
hereafter BCP), we explained how Fe can be present 
above the accreting NS photosphere. The larger entropy 



of the photosphere forbids convective mixing from the 
deeper burning regions (Joss 1977), making it unlikely 
that Fe is dredged up. In addition, the strong surface 
gravity rapidly depletes Fe from the upper atmosphere 
via sedimentation 1 . However, continual accretion during 
the burst at fairly modest rates (M > 10~ 12 M Q yr~ 1 ) 
can overcome Fe sedimentation. For EXO 0748-676, the 
persistent flux implies M > 2 x 10 _10 M Q yr^ 1 (Gottwald 
et al. 1986; CPM) and the sub-Eddington burst luminos- 
ity is unlikely to halt the accretion flow (see observations 
of GX 17+2; Kuulkers et al. 2002). 

If accretion is the source of the Fe that generates the 
absorption line, then accretion must be spread over a 
large fraction of the star. Since the NS surface is in- 
side the last stable orbit for the implied radius of R = 
AAGM/c 2 (CPM), the protons may reach the NS surface 
with a kinetic energy of w 200 — 300 MeV/nucleon, and 
decelerate via Coulomb scattering (Zel'dovich & Shakura 
1969). This filters ions by their charge and mass, stop- 
ping incident Fe above the photosphere at a column 
of « lgcm~ 2 (Bildsten, Salpeter & Wasserman 1992). 
The Fe remains there until it is either destroyed by in- 
cident protons (for M > 10~ 13 M Q yr -1 ) or sediments 
out. The deposition and nuclear destruction of Fe both 

1 Radiative levitation is a possible mechanism by w hich the Fe 
is suspended (BCP), an issue we will address in § 13.31 
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depend linearly on M and hence the steady state Fe 
column is independent of M when nuclear destruction 
dominates. For accretion of solar metallicity material, 
N Fc 3.4 x 10 19 cm~ 2 (BCP). While this Fe column is 
similar to a uniformly abundant solar metallicity atmo- 
sphere (Np c ~ 5 x 10 19 cm -2 ), the accretion/spallation 
scenario also predicts a plethora of elements with Z < 26 
produced by proton spallation of Fe (BCP, Table 1). 

In this paper, we discuss the physics of line trans- 
fer through a thin Fe layer above the continuum pho- 
tosphere. Previous theoretical work on the formation 
of spectral features focused on the Fe Lya line and 
edge (Foster et al. 1987; Day, Fabian & Ross 1992). 
However, these calculations presumed LTE and uniform 
abundance. In §|21we show that LTE is a poor assump- 
tion as the radiation field dominates the level popula- 
tion and ionization balance (London, Taam & Howard 
1986; BCP). We also discuss the line broadening physics 
and show that Stark broadening dominates for Fe Ha. 
We calculate in § O the Fe column needed to produce 
the observed Ha line and compare it to the Fe column 
from accretion/spallation. We also estimate the equiva- 
lent width of the resulting Lya and Pa transition, and 
constrain the rotation rate of the NS from the observed 
linewidth. We close by estimating the Fe column that 
might be sustained by radiative levitation, and summa- 
rize our results in §0] 

2. MODEL ATMOSPHERES 

We model the thin Fe line forming region as a one-zone 
scattering layer (Mihalas 1978) at fixed temperature and 
density that is illuminated from below by the contin- 
uum photosphere (see Figure We consider radiation 
and gas temperatures of fcsT g w fc^Tcff w 1 — 2 keV, 
in which case the atmospheric microphysics is simplified. 
The continuum opacity, /c con t = ^abs + ^Th, where K a b s 
is the absorption opacity (mostly free-free) and KTh is 
the Thomson scattering opacity, is dominated by elec- 
tron scattering. Bound-free transitions are suppressed 
since most atoms are fully ionized (Pavlov, Shibanov & 
Zavlin 1991; London et al. 1986). We do not consider 
the global problem of the thermal spectra from a burst- 
ing NS (for a review see Lewin, Van Paradijs & Taam 
1993 and references therein), but rather the simplified 
problem of line transport in the Fe layer. The spectra 
from a bursting neutron star, whose opacity is scattering 
dominated, is expected to be non-Planckian (Zeldovich 
& Shakura 1969; Castor 1974; Madej 1974; van Paradijs 
1982). For our purposes, we take the color correction 
given by Madej, Joss & Rozanska (2004) (also see Madej 
1991) to relate the observed color temperature to the ef- 
fective temperature. The effect of the color correction is 
nearly offset by the measured redshift correction. Hence 
the observed color temperature at infinity is roughly the 
effective temperature at the surface. For the sake of sim- 
plicity, we approximate the continuum flux with a black- 
body at the effective temperature. However, we will dis- 
cuss the effect of the spectral shape of the continuum flux 
on our results in § 13.11 

We assume that the continuum is formed below the line 
forming region and is freely streaming. The equation of 
radiative transfer is then (Mihalas 1978) 

/"-K— = ~ («abs + «Th + 4>l4>u) Iv + Kab S Bv + K T hJ V 
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Fig. 1. — Schematic model of the line forming region in a NS 
atmosphere during a Type I burst. 



+eip l 4> u B l/ + (l-e)ipi J R{v, v')J v ,dv\ (1) 

where I v is the frequency and angle dependent specific 
intensity, J v = (1/2) J I u dfi is the angle-average mean 
intensity, B u is the Planck function, tpi is the line opac- 
ity properly normalized, e = K a bs/^tot is the absorbed 
fraction, where K to t = K CO nt + V'ii 4> ls the line profile dis- 
tribution and is determined by the physics of § 12.11 and 
i?(^, v') is the redistribution function. The coordinate 
y measures the column from the base of the scattering 
layer. 

2.1. Line Broadening Physics 

We begin by discussing the line broadening physics. 
The number density of fully ionized H at the photosphere 
is 

where g is the surface gravity. It is sufficiently dense 
so that the intrinsic broadening of the hydrogenic Fe 
Ha line is dominated by the linear Stark effect (Paerels 
1997). The linear Stark effect is only relevant for hydro- 
genic atoms and breaks down for multi-electron atoms 
where the degeneracies of opposite parity states are bro- 
ken (Bethe & Salpeter 1957). The distribution of the 
electric microfields is given by the Holtsmark distri- 
bution (Mihalas 1978; also see Potekhin, Chabrier, & 
Gilles 2002) and the characteristic scale is given from in- 
teractions with neighboring protons, E = e/r^, where 
ro = (47rn p /3) -1 / 3 = 5 x 10~ 9 p^ 1 / 3 cm is the mean ion 
spacing and p = n p m p is the density. Hence the scale of 
Stark broadening is 

h 2 n 2 

A£stark = "Al/stark ~ Ty 9" 

m e Zrg 



Neutron Star Line Transport 
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for the Ha line, 0.6 eV for the Lya line, and 2.7 eV for the 
Pa line, where Z26 = Z/2Q. Note that we have chosen 
the upper state because of the strong dependence on n. 
The scale of thermal Doppler broadening is 



AE th = ftAi/th 
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where -Etran is the line energy, scaled to the Ha transition 
and A 56 = A/56. For Lya and Pa, AE th w 1.4 eV and 
A-Bth ~ 0.087 eV respectively. Stark broadening domi- 
nates the Ha line and Pa, but thermal Doppler broad- 
ening dominates the Lya line. The relative strength of 
Stark and thermal Doppler broadening scales like 



A = 



tark 



11 Z 26 A 26 



1/2 



IkeV 
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VlO 23 cm- 3 / 



-2/3 \n 



(5) 



where n is the upper state and n' is the lower state. Note 
the dependence on n and n' . Depending on the values 
of n and nf, the relative strengths of the two broaden- 
ing mechanisms can vary tremendously. For instance for 
the Ha transition, A ss 0.17, while for the Lya and Pa 
transitions A«2 and « 0.033 respectively. 

The timescale of changes in electric microfields is 
Tei = r /vth, where v t h — y / 2kBT g /m p is the ther- 
mal speed of the perturbing protons. At photo- 
spheric temperatures and densities, this gives tei 



7 x 10^ 16 s(fc B r 9 /lkeV)~ 1/2 (n p /10 23 cm- 3 )~ 1 / 3 , much 
smaller than the n — 3 — > 2 radiative deexcitation 
timescale of 3.1 x 10~ 14 s. Since the electric held changes 
during the time in which the atom absorbs and re-emits 
a photon, the energy of the two Ha photons are uncor- 
rected. With this in mind, we can assume complete re- 
distribution, greatly simplifying the line transport calcu- 
lation. 

The final broadening profile is a convolution of Stark 
and thermal Doppler broadening. In Figure [21 we 
show these profiles for various values of A. Asymp- 
totically, the distribution approaches the power law oc 
(A^/Az/stark) -2 ' 5 , expected for the Holtsmark distribu- 
tion (Mihalas 1978). For A 3> 1, we recover the thermal 
Doppler profile. For A <C 1, we approach the Holtsmark 
distribution. Note that for small values of A, there is a 
severe deficit at zero shift because this shift corresponds 
to perfect symmetry in the proton distribution. In fact 
for A = 0, the probability of a zero shift goes precisely to 
zero. The effect of thermal Doppler broadening is to fill 
in the distribution near zero (see Figure The deficit 
near zero shift yields an emission-like feature there, which 
must be interpreted as reduced scattering. This is simi- 
lar to the "ghost of Lya" where reduced scattering in the 
line contributes to an emission-like feature in the spectra 
of broad absorption line (BAL) quasars (Arav 1996). 

2.2. NLTE Effects 
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Fig. 2. — Distribution of frequency shifts due to a combination 
of thermal Doppler (Voigt) and Stark shifts (Holtsmark) at various 
relative strengths of each. The x-axis is given in units of AiAjtark 
and the relative strengths of each distribution is given by A = 
^Stark- The limits of A — and A 1 gives the Holtsmark 
distribution and the Voigt profile respectively. 



We now show that the atomic level populations in the 
photospheres of bursting neutron stars are not expected 
to be in LTE (London et al. 1986, BCP). For simplicity 
we only consider a two-level system. The rate equations 
for the density of upper (n u ) and lower (ni) states are 
then, 

= ni [Ciu + B lu J(E lu )} - n u [C ui + A ul + B ul J(Ei u )} , 

where Ci UlU i are the collisional rates from the lower to 
upper state and vice versa, A and B are the Einstein 
coefficients, Ei u is the energy of the transition and J is 
the mean radiation intensity at that energy. We take the 
dipole approximation from Jefferies (1962) for optically 
allowed transitions for the collisional rates: 



Ciu = 5.5 x 10 1 



f Ei u 
\k B T g 
-1.5 



exp 



Eiu 



IkeV 



. 10 23 < 



ksTg 

3) ft« S ' 1 



(7) 



where fi u is the oscillator strength for the transition. 
Putting this all together, we have C/„ ~ 2.2 x 10 9 s _1 
for the n = 1 — > 2 transition and 1.3 x 10 12 s _1 for the 
n = 2 — > 3 transition at ksTg = 1.8 keV. If we compare 
this to the radiative excitation/deexcitation rates in the 
same conditions, we find that Bi u J(Ei u ) = 7 x 10 13 s _1 
for the n — 1 — * 2 transition and Ri u = 7 x 10 12 s _1 , eas- 
ily exceeding the collisional rates for the n — 2 — > 3 tran- 
sition at a radiation temperature of ksTcg = 1.8 keV. 

Having radiative rates exceed collisional rates is not the 
only prerequisite to break LTE. If the radiation field were 
isotropic and Planckian, we would recover the Boltzmann 
distribution for the level populations where the relevant 
temperature setting the level populations is that of the 
radiation. Hence to break LTE, we require the radia- 
tion field to be either non-Planckian, anisotropic or both. 
This is indeed the case above the photosphere where the 
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radiation field is free-streaming and the source function 
deviates from Planckian due to the nature of scattering 
dominated atmospheres (Mihalas 1978). Therefore, to 
successfully model the spectral line formation, we must 
solve the equations of radiative transfer and statistical 
equilibrium simultaneously. 

3. LINE FORMATION CALCULATIONS 

We rewrite equation by defining the flux F v = 
(1/2) f pl„dp, and K v = (1/2) J p?l v dp. Integrating 
equation (JTJ over dp and pdp we find (Mihalas 1978), 

OF 

v =-j u + {i + [(i-p+ ePv) b v + P j„ 



+ (l-e)(3 v R{v,v')J v ,dv' 



dK v 



= -F v , 



(8) 
(9) 



where dr v = (K abs + n Th + ipi<p„) dy, p = 

«Th/ («abs + KTh) and (5 V = ^l4> v / (tabs + KTh). 

Note that t„ is defined from the base of the scattering 
layer. If we define the variable Eddington factor as 
/„ = K v /J v , we can write equation QjJJ as 

<> ' J '" T " Ju~(l + PvT 1 [0--P + <&0 Bv + P J V 



+ {l-e)j3 v / R{v,v')J v ,dv' 



(.10) 



Equation I|1U|) is a generalization of Harrington's (1973) 
equation including the use of variable Eddington factors. 

We solve equation (|10[1 with a standard methodology 
(Mihalas 1978). Following the nomenclature of Mihalas 
(1978) wc define a symmetric and antisymmetric average 
for the intensity 

Vv{TuM) = n (11) 



v v (t v ,\(i\) = -\I v {t v ,\h\) - I v {t v ,-\h\ 



(12) 



which combined with the equation of radiative transfer 
(eq.pQ) yields 

\p\-^— = -U v + b v , 



<9t„ 
,du v 



where 



S v = (1 + PJ) 1 [(1 - p + epv) B v + P J V + 
(l-e)j3„ / R{v,v')J v ,dv' 



(13) 
(14) 

(15) 



is the source function. To solve this system of equations, 
we take an initial guess for the variable Eddington factor, 
/„, and assume an initial input I v at the bottom of the 
scattering layer that is isotropic in the outward half space 
(p > 0). Wc then solve for J„ in equation (|l()(l . We 
compute S v given J v and solve for u v and v v in equation 
H13(l . With u v and v v in hand, we compute, /„, F v , and 
K Vl compute a new f v and iterate to convergence. 

Since the atmosphere is in NLTE and scattering dom- 
inated, we must simultaneously solve the equations 



of radiative transfer (eq.pQ) and statistical equilibrium 
(eq.jHJ). For simplicity, we assume a two- level atom and 
calculate the level populations between the two states 
(in the case of Ha, the two states are the n=2 and n=3 
states; for Lya, n=l and n=2; for Pa, n=3 and n=4) 
fixing the sum of the two states. We iteratively solve 
the radiative transfer equation (eq.pP) and use the com- 
puted mean intensities to solve for statistical equilibrium 
(eq. |S] ); i.e. A iteration. We define convergence when the 
mean intensity varied by less than 8 J/ J < 10 -3 , which 
was typically achieved in 3-5 iterations. 

3.1. Line Profile and Equivalent Width 

We solve the thin scattering layer illustrated by Figure 
^with the above technique. However, before we continue, 
we review the analytic result of a specific case of scat- 
tering dominated atmosphere known as Schuster's law 
(Schuster 1905). Assuming that line opacity is dominant 
(fi v 1), the absorbed fraction is small (e — > 0), and 
the scattering is coherent {R{v\v) = 5 U ' rU ), equation 
ifTHjl becomes 

d 2 f J 

W = a (16) 

Solving this equation, we find f v J v — A v t v + B v , where 
A v and B v are constants. From equation JJjJ, wc also 
find F v =f v A v . 

We now make some further approximations to deter- 
mine A v and B v . We assume the Eddington approxi- 
mation everywhere so that f v = 1/3 (Rybicki & Light- 
man 1979). We also assume the two stream approxi- 
mation, which breaks the radiation field into one outgo- 
ing stream = I v (t v , \p\)) and one incoming stream 
(/+ = 1 v (t Vi — \p\f). Consistency with the Eddington ap- 
proximation demands \p\ = 1/V3 (Rybicki & Lightman 
1979). The boundary conditions are that the incoming 
flux at the surface is zero, I v (t v — r* ot , p = — 1/\/3) = 0, 
where r* ot is the total optical depth of the scattering 
layer at frequency v and that below the line forming re- 
gion, the incident intensity is fixed at some value typ- 
ical for Type I bursts {I u {t v = 0, p = l/\/3) = iinc)- 
Per frequency the optical depth of the slab is given by 

r*°V(V), where r*° 
optical depth and is given by 



7"g 0t — iVFc,n=20'o,23 

iVpie n =2 
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where ero,23 = 7!"e 2 /23/(?Ti e cAi'stark) is the cross section 
at line center. In the two-stream approximation, our 
definitions for the mean intensity and the flux become 



F u = 



2V3 



and hence solving for A v and B v , we find 
j _ - 2 

-*■ trans — 



V3t 



(18) 
(19) 

(20) 



where /trans = I v {t v = T* ot ,^i = l/\/3) is the trans- 
mitted intensity. Equation l|20|l is known as Schuster's 
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Fig. 3. — Hydrogenic Fe Ha line profiles of various values of e in 
the line at the NS surface. For the background we take Tq * = 100, 
^B^eff = 1.8 keV and n p = 10 23 cm -3 . The numerical cal cula tion 
of two stream approximation and exact analytic result f eg 1201 ) arc 
indistinguishable and are represented by the thin-solid line. The 
full numerical results properly angle averaged for complete redis- 
tribution are represented by the other curves for various values 
of the absorption ratio e. The difference between the thin-solid 
line and the thick-solid line are due to the departures of the two- 
stream approximation compared to the full angle-averaged solution 
and coherent scattering compared to complete redistribution. The 
difference between the other curves are due to different values of e. 



law. Note that at large T* ot , the specific intensity scales 
with l/r* ot unlike the absorption dominated case where 
I v — ► B v at large r* ot . This difference in scaling with op- 
tical depth between a scattering dominated atmosphere 
and an absorption dominated atmosphere is the key to 
understanding the nature of the Fe Ha line. 

We compare the flux from our exact numerical cal- 
culations with Schuster's law (eq.[5Dj) in Figure |3] For 
simplicity of the discussion, we take a simple model of 
a hydrogenic atom whose energy levels are defined by 
E n = — 13.6 Z 2 n~ 2 eV, which gives the energy of the Fe 
Ha transition of 1.2769 keV. We will consider the fully 
relativistic atom with the complications of fine structure 
splitting later. We take rg* = 100, k B T cS = 1.8 keV 
and n p = 10 23 cm -3 for these models. Our numerical 
calculation for coherent scattering in the two-stream ap- 
proximation and Schuster's law are indistinguishable and 
are both represented by the thin-solid line. This con- 
firms the accuracy of our numerical treatment with the 
exact analytic result. The other calculations represent 
the full numerical solutions with complete redistribution 
with the variable Eddington factor, /„, computed self- 
consistently. These curves are computed with 10 angular 
points, 100 frequency points and 100 depth points at var- 
ious values of the absorption coefficient e. The difference 
between the thick-solid line and the thin-solid line is due 
to differences between coherent scattering and complete 
redistribution and the nature of the two-stream approx- 
imation which fixes /„ = 1/3 at every point compared 
to the full self-consistent calculation. Note the emission- 
like feature at line center, which is due to the deficit near 
zero shift discussed in § 12.11 



Now that we have established the effectiveness of our 
numerical code, we discuss the nature of the hydrogenic 
Fe Lya, Ha and Pa lines. Since we have limited our 
study to two-level atomic systems, the relevant transi- 
tions starting from the Is ground state are Is <=^ 2p, 2p 
3d and 3d 4f. The fine structure splitting of these 
transitions significantly complicates the line structure 2 . 
Due to fine structure splitting the energy levels of an 
hydrogenic atom are 



Z 2 

E n<j = -13.6^ 
n 2 



1 



(aZ) 2 



_3_ 

An 



eV, (21) 



where j is the total angular momentum quantum num- 
ber, 



k = 



I 



if 3=1 



l + l iij = l + 



(22) 



I is the orbital quantum number, and a is the fine struc- 
ture constant. The Lya Is 2p transition is fine struc- 
ture split into two lines with relative line intensities of 2:1 
with energies of 6.9728 keV and 6.9521 keV, a difference 
of 20.7 eV. The Ha 2p <=^ 3d transition is split into three 
lines with relative line intensities of 9:1:5 and energies 
of 1.2811 keV, 1.2790 keV, and 1.2997 keV. Taking the 
energy of the strongest line to be the centroid, the split- 
ting of the levels are 0, -2.1, 18.6 eV (Bethe & Salpeter 
1957). The size of the splitting is comparable to the EW 
and the FWHM of the observed line. Finally the Pa 3d 
4f transition is split into three lines with intensities 
20:1:14 and energy differences of 0, -0.43, 1.61 eV. These 
splittings are smaller than the effect of Stark broadening, 
so that their effect on the line profile is small. 

We illustrate the line profile with fine structure split- 
ting in Figure 0] We also turn on statistical equilibrium 
so that from an initial state where the initial popula- 
tion of the two-level atom was assumed to be in the n=2 
state (with an optical depth normalization of Tq ' = 100), 
we simultaneously solve radiative transfer and statisti- 
cal equilibrium. Since the opacity is dominated by line 
itself, the absorbed fraction is determined from the ra- 
tio of the collisional rate to the total transition rates 
(e = Kabs/^tot ~ Ci u /(Ri u + Ciu)). In statistical equi- 
librium, e relaxes to a small value of 0.033, typical of a 
scattering dominated atmosphere. This line profile rep- 
resents the complete solution to the spectra from the 
resonant scattering layer. Note that the triplet structure 
of the fine structure lines is reduced to a doublet because 
Stark broadening is so large. 

The line's equivalent width (EW) is the proportional 
deficit and is defined as 



EW 



1 



F, 



line 



^cont. 



du. 



(23) 



Though rotation will modify the line profile, it con- 
serves equivalent width (Gray 1992) to first order in 
QR/c, where f2 is the angular rotation rate. In Figure 
|SJ we show the fractional equivalent width, defined as 
the equivalent width in wavelength units over the line 
center wavelength, as a function of the Fe column in 
the n — 2 state for n p = 5.2 x 10 22 cm~ 3 (dotted line) 
and n p = 10 23 cm~ 3 (solid line), corresponding to the 

2 We thank Rashid Sunyaev for raising this important issue. 
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Fig. 4. — Hydrogenic Fe Ha line profiles with fine structure split- 
ting. We take the same background as Figure 131 ( k nT„g = 1.8keV 
and n p = 10 23 cm - 3 ). We also simultaneously solve statistical 
equilibrium assuming an initial two-level system with everything 
in the n=2 state with an optical depth normalization of Tg ot = 100. 
With these parameters, statistical equilibrium gives e 0.033. 
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Fig. 5. — Fractional Equivalent width as a function of the Fe 
column in the n = 2 state for two different proton densities of 
n p = 5.2 X 10 22 cm -3 (dotted line) and n p = 10 3 cm -3 (solid line) 
at fcsT c fj = 1.8 keV. We show the observed EW as a banded region 
and the intersection of these lines with the banded region indicates 
the n = 2 Fe column needed. We also present an extremal LTE 
calculation (dashed line) where Aq Rj 0.2 (BCP) for comparison. 



Coulomb stopping density for radial infall and the r w 1 
surface (approximating a solar metallicity atmosphere) 
for a radiation temperature of ksT c s — 1.8 keV. We 
also show the observed EW (CPM). The scattering dom- 
inated calculations imply A/p e ,n=2 ~ 5 x 10 17 cm~ 2 . We 
also plot an extremal grey atmosphere model in perfect 
LTE (e = 1) for comparison. The LTE calculation re- 
quires a much larger N-p e n=2 and highlights how poor 
the LTE calculation fares compared to the full resonant 
calculation. 
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Fig. 6. — iVF e ,n=l) N-p e , EWLya, and EWp as a function of 
total EWhq. for parameters n v = 5.2 X 10 22 cm~ 3 (thin lines), 
n v = 10 23 cm -3 (thick lines) and fcsT ef j = 1.8 keV. The banded 
region represents the observed EW and the intersection of these 
lines with the banded region indicates the Fe column (upper plot) 
needed in the hydrogenic ground state (dashed lines) and total Fe 
(solid lines) and the predicted EW in the Lya line (solid lines; 
lower plot) and Pa line (dotted lines; lower plot). The Pa values 
have been divided by 10 so that they may fit on the same plot. 



We calculate a grid of models for the Lya, Ha and Pa 
transitions to relate the equivalent widths and level pop- 
ulations in each transition to one another. We estimate 
the total ATpe to that in the n=l state via Saha equilib- 
rium. Since the radiation field drives ionization, we ex- 
pect our results from Saha equilibrium will give a lower 
bound on N Fe . We relate EW Hct , -/V Fe ,n=i, N Fe , EW LyQ 
and EWp Q in Figure [S] at two different proton densities: 
n p — 5.2 x 10 22 cm -3 (Fe stopping layer; dotted line) and 
n.p = 10 23 cm -3 (solid line). The larger proton density is 
indicative of what would happen in a uniform metallicity 
atmosphere. The observed Fe Ha EW gives an implied 
column of A^F e ,n=i « 1 — 2 x 10 19 cm -2 if all the Fe is 
concentrated in the Fe stopping layer. Saha equilibrium 
then implies a total Fe column of Np c = 3 — 6 x 10 20 cm~ 2 
at the Fe stopping layer, a factor of « 10 higher than the 
solar metallicity accretion/spallation scenario. For the 
larger proton density (n p = 10 23 cm -3 ), the implied Fe 
column is w 10 20 cm~ 2 , a factor of two larger than that 
implied by a uniform solar metallicity photosphere. Since 
ionization balance is likely given by the radiation field 
(London et al. 1986), the true Np e is likely even higher 
for both these scenarios. We also estimate that our im- 
plied Nfb,b.=i gives an EWL yQ » 15- 20 eV (redshifted 
11-15 eV). Though this is larger than EWh q , rotational 
broadening has a large impact on its detectability (see 
§E2J. For the Pa line, the EW Pa » 4- 7 eV (redshifted 
3-5 eV) for the observed EWh q - 

We now take ksT e ff as a free parameter to under- 
stand how our results depend on the color correction 
and spectra. We plot the implied n=l and total N Fe 
as a function of ksT e s at the observed EW in Fig- 
ure \7\ Since we are interested in Fe at the stopping 
column (measured from the top of the atmosphere) of 
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Fig. 7. — The Fe column required to generate the observed EW 
in the Fe Ha line as a function of fcsT e ff at the Fe stopping column 
of 2/stop ~ lgcm -2 (thin lines) and jfetop = 2 gem - 2 (thick lines), 
which is more representative of a uniform atmosphere. The various 
Fe columns represented are the hydrogenic Fe column (dottcd-linc) , 
fully ionized Fe column (dashed-line) and total Fe column (solid- 
line). We find that TVpe must at least be a factor of 3-4 larger 
than that predicted from the accretion/spallation scenario of 3.4 X 
10 19 cm - 2 (BCP), but is within a factor of 2 of a uniform solar 
metallicity atmosphere of 5 X 10 19 cm -2 . 
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Fig. 8. — Redshifted line profiles for the Fe Pa, Ha, and Lya 
at various rotation rates. We fixed the Fe Ha EW to the ob- 
served values and assume that the Fe is concentrated in the accre- 
tion/spallation stopping layer. We compare models (from bottom 
to top) at u sp i n = (1 + z) -1 Sl/(27r) = Hz (thick-solid lines), 
45 Hz (short-dashed lines), 100 Hz (dotted lines), 200 Hz (long- 
dashed lines), and 300 Hz (dashed-dotted lines) for sini = 1. We 
do not include a 300 Hz Lya line profile as it is totally rotationally 
broadened away. 



J/stop ~ 1 g cm 2 (thin lines) , the local proton den- 
sity, which is important for Stark broadening, is then 
n p » 9.4 x 10 22 (fcsTcff/lkcVr 1 cm" 3 . The implied Fe 
column is minimized at ksT e s = 1.3 keV, but it is 3- 
4 times larger than expected from our solar metallicity 
accretion/spallation scenario. We also plot a similar Fe 
column at y s t op — 2 gem -2 (measured from the top of 
the atmosphere), which is more representative of a uni- 
form atmosphere. From Figure HO and we find that our 
the implied Fe column for a solar metallicity atmosphere 
is within a factor of 3 for the entire range of burst tem- 
peratures of 1-2 keV. We discuss these points further in 

3.2. Rotational Broadening 

Rotational broadening heavily modifies the structure 
of the lines. We estimate its effects on the line profile as 
follows. Since the photon occupation number n oc I v /v 5 
is invariant under Lorentz transformations, the effect of 
rotational broadening on the emergent NS flux is 



jpobs 



1 

27 



271 



d(cos 9) cos 9 



(1 — (ilR/c) sin 9 cos <f> sin i) 

I (v(l — (QR/c) sin 9 cos <j> sini), 

cos 9(1 + (flR/c) sin 9 cos 4> sini)) . 



(24) 



Following various authors (Chandrasekhar 1945), we ig- 
nore the effect of Lorentz boosts on the angle and focus 
on the effects of the frequency shift. We are currently 
calculating the exact general relativistic calculation else- 
where (Morsink et al. 2005; also see Ozel & Psaltis 2003 
and Bhattacharyya, Miller & Lamb 2004). In addition 



for narrow lines, \(dl v / d\\iv)(Q,R/ c)\ ^> 1, so we can ig- 
nore the effects of the (1 — HR/c) 3 in the denominator, 
leaving 



pok 
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2tt pi 

d<j) / d(cos 9) cos ( 
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Jo 
ttR 



■ sm v cos sin % 



cos 9 , (25) 



showing that rotational line broadening goes like 
(flR/c) sini. We plot the Fe Pa, Ha and Lya lines for a 
1.4Mq NS with R=9. 2km at various rotation rates up to 
300 Hz in Figure|H] To ease comparison with observation, 
we plot the lines as a function of wavelength and include 
the redshift correction of z — 0.35. We also include the 
1 + z redshift correction to the rotation rate. For the 
highest rotation rates (> 200 Hz), the Fe Ha central line 
depth would only be about 5-10% below the continuum 
and hence the line would be difficult to detect at the ob- 
served EW. At lower rotation rates the doublet nature of 
the line gives an asymmetric profile with a sharp feature. 
For the Pa line, the fine structure splitting in the line is 
smaller than the Stark effect so it appears as a singlet. 

Rotational broadening is more pronounced on the Lya 
transition, which has a FWHM w 130 eV at a spin fre- 
quency of 45 Hz. This is 5 times larger than the rota- 
tional FWHM for the Ha line due to the larger energy 
of the Lya transition compared to the Ha transition. At 
the implied EWL ya matched to the observed EWh q , the 
line is extremely shallow with a central line depth of ap- 
proximately 10%. Hence even at this modest rotation 
rate, the Lya line may be difficult to detect. For larger 
rotation rates (> 100 Hz), Lya is so distributed that it 
may be impossible to detect. 

Observations of dips and eclipses on EXO 0748-676 
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(Parmar et. al. 1986) have revealed that the orbit of 
the secondary is nearly edge-on (sini w 1). Assuming 
the rotation axis of the NS is aligned with the orbit, we 
favor very low rotation rates (50—100 Hz) for EXO 0748- 
676 compared to the rotation rates that are typically 
observed in other LMXB systems where the rotation rate 
is typically 300 — 600 Hz (for a review see Strohmayer 
& Bildsten 2003). The recent observation of a 44.7 Hz 
burst oscillation (Villarreal & Strohmayer 2004) is thus 
consistent with the observation of the line, assuming that 
this oscillation is indicative of the NS rotation rate. It 
raises the exciting prospect of fitting the line with an 
appropriate spectral calculation to derive the NS radius 
(Morsink et al. 2005). Other accreting NSs may also 
possess atomic lines, but we would not expect them to 
be detected unless they rotate slowly (< 100 Hz) or are 
highly inclined. 

3.3. Radiative Levitation 

BCP noted that Fe in the atmosphere may be sus- 
tained by radiative levitation similar to what occurs in 
hot stars (Michaud 1970; Michaud et al. 1976) and white 
dwarfs (Vauclair, Vauclair & Greenstein 1979; Vennes et 
al. 1988; Chayer, Fontaine & Wesemael 1995). To be- 
gin, we calculate the radiative force on a hydrogenic Fe 
atom, Fj-ad = c _1 J a v F v dv„ from reflection of the inci- 
dent flux, F v , where a v is the cross section and is defined 
as a v = CTO4201/, Co, 12 = ve 2 /i 2 /(m e cAi/), <j) v is the con- 
voluted Holtsmark profile, and fn = 0.42 is the oscillator 
strength for the Lya transition (Bethe & Salpeter 1957). 
Assuming that the flux does not change as we integrate 
across the line, J a v F v dv = F„ire 2 j 12/ (m e c). For a ra- 
diation field that is isotropic in the outward hemisphere, 
F v = nB v (T) (Rybicki k Lightman 1979), B V {T) is the 
Planck function which we approximate with the Wien 
tail, we find. 



F r 



ad 



9 x 10 3 exp 



6.9 keV 



^56 9l4, J 



(26) 



Am p g 

where g u — g/10 14 cms~ 2 . For fcsT c ff = 1 keV, the 
radiation force on the Lya exceeds gravity by a factor of 
9 if all the Fe is hydrogenic. 

Since the radiation force exceeds that of gravity, the 
Fe will not necessarily sediment. Rather, Fe can be sus- 
tained in the atmosphere, allowing for a potential con- 
centration. As the column of Fe builds, more and more 
of the line flux is absorbed which decreases the radia- 
tive force. This continues until an equilibrium column is 
reached where enough of the line flux is absorbed that 
the Fe feels the same net force as a background pro- 
ton, so that both species have the same scale height 
(Chayer, Fontaine & Wesemael 1995). This is achieved 
when the radiative force is equal to the effective force 
-Frad = Am p g c ff = [A — (Z + l)/2] m p g, which takes into 
account the effects of both gravity and the electric field 
for a hydrogen atmosphere and assumes all the Fe is hy- 
drogenic (Vennes et al. 1988; Chayer, Fontaine & Wese- 
mael 1995). 

We now estimate the equilibrium column sitting above 
the continuum photosphere. The radiative force on a 
scattering layer is 



Frad = 

lines 



tEWi 



(27) 



where i/i is the energy of the ith line, -Fl^cont is the contin- 
uum flux and EWi is the equivalent width of the ith line. 
For pure coherent scattering, Schuster's law (eq.[20j) and 
the definition of the equivalent width (eq. [231) gives 



EWi = 



'OA ' 



2/V3- 



-dv, 



(28) 



'OA • 



where 



and 



is the normalization of the optical 



depth and line profile distribution of the ith line respec- 
tively. We find an extremely simple form for equation 
l|28|l by taking the asymptotic form 3 of the Holtsmark 



distribution tf> v w 0.75 {\v — ^o|/Ai/stark) 2 5 - Assuming 
T0 Ot is large, we break up the integral in equation l|28|l 
into a regime where Tq *^ is large compared to 2/\/3 
and a regime where Tq '^ is small compared to 2/-\/3. 
Thus, we get 



EW l ^2r^Av stalk 

«2.8Al* tark (T$) 



drj 

Itot 



0.757? 



-2.5 



'/o 



2/V3 



-dr\ 



(29) 



where t\ = \v- fo|/Az/ stark and 0.75^°% 2 5 = 2/>/3. 
We extend this simple form for the EW to handle fine 
structure splitting if the split lines are widely separated 
and therefore distinct. We sum over each split line to get 



EWi » 2.8Az/ St ark (t3 



,0.4 \ - 



(30) 



where Sij are the relative strength of each split 
line for the ith transition. For Lya, which is 
split into two lines with a ratio of 2:1, we find 

EW Lya = 4.2AiA3tark,Lya (^Lya) ' 4 ■ Similarly for 
Ha with a line strength ratio of 9:1:5, EWn a = 
5.2Ai/ Stark , HQ (r*°t Q ) a4 . 

Following our previous discussion, we now equate the 
radiative force fea.|27|) to the effective force to find the 
equilibrium column. For simplicity, let us first consider 
only the Lya line and a radiation field that is isotropic 
in the outer half hemisphere. We find 



^Fc.tot « 5 x 10 2 



3 x 10 14 cms 



5/3 



" 2 / 3 cm- 2 , 



(31) 

for kT 2 keV, where £ = iVF e ,tot/A?Fe,n=i is the ratio of 
total Fe column to ground state hydrogenic Fe column. 
The Lya equivalent width that is associated with this Fe 
column is 



EW^ya ~ 400 



9 



3 x 10 1 



2/3 



C 2/3 eV. (32) 



Radiative levitation can easily sustain an atmosphere 
whose metallicity can be supersolar. In Saha equilibrium, 
the ratio between fully ionized Fe and the n=l hydro- 
genic state is 10-20 and A FC:tot = 7 - 12 x 10 19 cm~ 3 , 
larger than a constant solar metallicity atmosphere 

3 The standard definition is <j> v 1.5 {\u — wo\/ A^g tar k)~ 2 ' J for 
large (\v — i^o|/Afg t ark)- The numerical prefactor of 0.75 comes 
from the fact that the line is symmetric and hence the distribution 
runs from —00 to 00 and not from to 00 
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TVpc ~ 5 x 10 19 cm~ 2 (BCP). These values are also tanta- 
lizingly close to what would be required to for the mea- 
sured Ha line, but we note that the radiatively levitated 
Fe column is sensitively dependent on the background 
continuum and surface gravity. Since bursting NS atmo- 
spheres are expected to be harder than a blackbody as 
discussed earlier in § |3 we would expect these values to 
be lower limits. The associated Lya EW for these Fe 
columns in this approximation is 50-90 eV. 

Other transitions also contribute to the radiative force 
on the Fe column. We now show that their contributions 
are small compared to the Lya radiative force. First, we 
consider the contribution from the Ha transition. For a 
Boltzmann distribution between the n=l and the n=2 
state at kT — 2keV, the column density in the n=2 
state is -/Vp Cjn= 2 = 0.1AF C . n =i- The radiative force from 
the Ha transition is -F r ad,Hc« ~ 2 x 10 12 g cm s~ 2 which is 
a factor of 5 smaller than the radiative force on the Lya 
transition of F r ad,Lya ~ 10 13 g cm s~ 2 . The contribution 
from the Ha transition increases the column by about a 
third, but we will ignore this contribution in keeping with 
the order of magnitude spirit of this calculation. Other 
transitions also contribute, but their radiative forces are 
weaker as their fluxes are weaker by the Rayleigh- Jeans 
factor of v 2 compared to Ha. 

4. DISCUSSION AND CONCLUSIONS 

Motivated by the accretion/spallation scenario that 
we first presented in BCP, we have calculated the res- 
onant radiative transfer of a thin scattering layer above 
the NS photosphere. While our calculation and micro- 
physics are geared toward understanding the feature as 
an Fe Ha line, the basic physics would remain the same 
for other atomic transitions. Namely, if the line is pro- 
duced near the continuum photosphere where the den- 
sities are relatively high, Stark broadening is dominant, 
NLTE effects are prevalent and the line is a result of 
resonant scattering. Applied to Fe Ha, our calculations 
show that the EW of the Fe Ha line observed by CPM 
from EXO 0748-676 requires an Fe column 3-10 times 
larger than the accretion/spallation scenario predicts for 
accretion of solar metallicity material (BCP). We also 
approximate a uniform solar metallicity atmosphere by 
assuming the Fe sits at a larger depth. In this case, we 
find that the required column of Fe is within a factor of 
3 of a solar metallicity atmosphere over a large range of 
temperatures. 

The fact that the Ha line is narrow requires that the 
rotation rate of EXO 0748-676 is slow. Larger spin rates, 
assuming a neutron star radius of around 10 km, would 
wash out such a line unless the surface emission is con- 
fined to a narrow region around the spin axis (Bhat- 
tacharyya, Miller & Lamb 2004). A slow rotation rate 
for EXO 0748-676 is consistent with the recently mea- 
sured 44.7 Hz burst oscillation (Villarreal & Strohmayer 
2004). With these line profiles in hand and assuming 
that the burst oscillation of 44.7 Hz is the rotation rate 
of the NS, detailed fits to the data raises the exciting 
prospect of fitting for the neutron star radius (Morsink 
et al. 2005). 

Since the fine structure splitting of Fe Ha is of order the 
size of the rotational broadening, we cannot assume that 
the line is narrow compared to the rotational broaden- 
ing. Previous work (Ozel & Psaltis 2003; Bhattacharyya, 



Miller & Lamb 2004) has presumed that the intrinsic line 
profile does not matter because rotational broadening is 
dominant. While true for NS rotating at 300-600 Hz, it 
is not the case for a Fe Ha line at 45 Hz for a 10 km NS. 
Rotational broadening may be the largest effect, but it 
is by no means the dominant effect. Therefore, the in- 
trinsic line profile must be known to apply the observed 
lines to the study of NS physics. 

We also estimated that the Fe EWL ya and EWp Q are 
approximately 15-20 eV and 3-7 eV respectively at the 
observed EWh - We can extend this calculation to in- 
clude other n=l and n=2 transitions aside from Lya and 
Ha. However, because of the greatly reduced oscillator 
strengths, these lines are significantly weaker. For in- 
stance, we have modeled H/3 EW to be < 2 eV. Rota- 
tional broadening makes a huge impact on the observ- 
ability of the Fe Lya transition. While the spectral res- 
olution of Chandra and Astro E2 is more than adequate 
to detect such a line, the Lya line will be shallow and 
difficult to detect. We expect rotational broadening to 
be a small effect for Pa, however, the redshifted energy 
of the Pa transition of 330 eV would make it difficult 
due to interstellar absorption. 

In calculating these profiles of Stark broadened lines, 
we have assumed a pure hydrogen atmosphere. However, 
the Stark broadening scale given in equation Q depends 
on the background composition. Namely Ai?g ta rk °c 
(2A/(Z + 1)) 2/3 Z 1 / 3 . For solar composition, A£ S tark 
is larger by 14%. This reduces the required N-p e needed 
to reproduced the observed EW by 10%. For a pure He 
atmosphere A£stark is larger by a factor of 2.42, which 
reduces the required Np e to reproduce the observed EW 
by 60%. For an atmosphere that consist of equal parts 
He and H by number, we find the required N-p e is reduced 
by about 50%. Such He-rich atmosphere would demand 
accretion from a H-poor donor such as 4U 1820-30, but 
the initial spallation would reduce the He fraction at the 
photosphere. 

The formation of the He-like Fe spectral feature in the 
low-temperature part of the burst is an issue that war- 
rants further study. In Saha equilibrium, H-like Fe tran- 
sitions to He-like Fe at 1.2 keV (CPM). We note however 
that the microphysics of this line is significantly differ- 
ent from the Ha line. Its instrinsic broadening is due to 
thermal Doppler effects because the lack of degenerate 
opposite parity states suppresses the linear Stark effect 
(Bcthc & Salpeter 1957). Modeling this feature requires 
a separate calculation which is beyond the scope of this 
paper. 

If accretion/spallation is indeed not the endpoint of 
accretion onto a neutron star, we expand on radiative 
levitation as a means to suspend the Fe, which we first 
pointed out in BCP. We give a simple estimate for the Fe 
column that could be suspended by radiative levitation. 
We find that radiative levitation could suspend a column 
of Fe needed to produce the line. However, this result de- 
pends sensitively on the exact nature of the background 
spectra. Further work will be necessary as the contin- 
uum spectra could be much harder than the blackbody 
background that we have assumed. 

Though we assume solar metallicity throughout our 
calculations, the composition of the accreting material 
remains highly uncertain. If the metallicity of the flow is 
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sufficiently supersolar, the disparity between the accre- 
tion/spallation Fc column and that implied from obser- 
vation can be resolved. However, the contrast between 
the deep dips and the unobscured state of EXO 0748-676 
indicates that the metallicity may be subsolar by a factor 
of 2-7 (Parmar et al. 1986). More recent work has found 
that the abundance ratios of Mg, Ne and O are consis- 
tent with solar with large uncertainties (Jimenez-Garate, 
Schulz & Marshall 2003). The intrinsic metallicity of the 
flow remains largely uncertain. 
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